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The entanglement in many-electron states is investigated using a global entanglement measure, 
viz. average site mixedness. We have examined metallic states of noninteracting electrons, Nagaoka 
and Gutzwiller states of strongly-correlated electrons, and superconducting states. Uncorrelated 
metallic states at half filling seem to maximize entanglement, as these states optimize the number 
of holes, the number of doubly-occupied sites. Entanglement is calculated explicitly for Gutzwiller- 
projected many-electron states in one dimension, which have less entanglement as double occupancy 
is inhibited in these states. Entanglement in superconducting states, which tend to promote double 
occupancy, is calculated as a function of the energy gap, and found to be lower than the metallic state 
entanglement. There is a possibility of a regime with a nonzero single-site concurrence depending 
on the energy gap. 



Quantum entanglement is perceived as a resource 
for quantum communication and information processing, 
and has emerged over the last few years as a major re- 
search area in various diverse fields such as physics, math- 
ematics, chemistry, electrical engineering and computer 
science[l|, 0, 0]. There has been a wide-spread investi- 
gation of entanglement properties of the ground states 
and excited states of well-studied quantum spin models 
of the condensed matter physics like the Heisenberg-XY 
models, the transverse-field Ising modeljl, H, B S B 111- 
Also, there has been an intense study of the dynamics of 
entanglement of these models, both from numerical and 
analytical approaches. 

The spin-only states have several advantages; many 
entanglement measures have been successfully used to 
characterize and quantify the entanglement of many- 
spin states of lattices. Pair-wise entanglement or 
concurrence [l^l, which measures the entanglement of a 
given pair of spins in a many-spin state, has been widely 
studied both numerically and analytically 0, 0, The 
global entanglement measure llj , which quantifies how 
entanglement is distributed and shared over various spa- 
tial parts of the lattice, has been studied. The reduced 
fidelity susceptibility measure has been studied to in- 
vestigate the critical properties of spin systems [l^]. In 
contrast, the entanglement properties of general many- 
electron states, which combine both spin and orbital de- 
grees of freedom, have not been studied extensively. In 
this article we will use a generalization of the global en- 
tanglement measure, viz. the average site mixedness, to 
quantify the entanglement between spatial parts of many- 
electron states on lattices. 

Let us consider electrons either localized or hopping 
around on TV lattice sites. The details of the interactions 
determine the exact many-electron ground state that we 
may like to study. Since the entanglement properties 
can be discussed entirely through the many-body state 
itself, we will refer to the interactions only through the 
structure of the state in this article. There are four states 
per site, implying a basis of 4^ states for the iV— site 
lattice. A convenient basis states for any site is given by 
|0), I T)i I i), I Ti); which correspond to no occupancy (or 



a hole), up-spin electron occupancy, down-spin electron 
occupancy, and double occupancy respectively. A general 
pure many-electron state for TV sites can be written as, 
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where labels the basis states of site i listed above. 
This state is characterized by 4^ complex-number am- 
plitudes {4'ai,a2..} for the basis states. The entangle- 
ment distribution in the state can be investigated us- 
ing various reduced density matrices (RDM) constructed 
through partial traces over the density matrix of the 
system p — \il]){ip\. In general a RDM of a subsystem 
would correspond to a mixed state, implying an entropy 
for the subsystem, and a bipartite entanglement between 
the subsystem and the rest of the system, and also some 
amount of multi-party entanglement between the vari- 
ous sites of the subsystem. It is quite an uphill task to 
characterize and quantify the entanglement distribution 
in various partitions of the system. 

In the context of many-electron states even two-site 
RMDs are difficult to work with and unravel the pair- wise 
entanglement, as the RDM is now 16-dimensional since 
no occupancy and double occupancy can occur. Here we 
will use the single-site RDMs, and use a generalization of 
the global entanglement measure [llj, that has been used 
extensively in spin systems. The entanglement measure 
is given in terms of the single-site RDM pi as, 
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In the above, the pre-factor is so chosen that if each 
site density matrix is maximally mixed (corresponding 
to all four eigenvalues being equal to 1/4), the entangle- 
ment measure equals unity. Since, 1- Trpf would be zero 
for a pure single-site RDM, the above measures average 
site mixedness or how entangled is a site with the rest 
of the system on an average. Though we are working 
with strictly local, single-site RDMs, the averaging pro- 
cess makes it a global measure. If a large number of sites 
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with nonzero measure are entangled with other sites, we 
get a nonzero entanglement. 

Let us consider a quantum state with — nN number 
of electrons, where n is the electron density. The state is 
further characterized by the densities of up and down spin 
electrons, with 7V| — n^N, Ni = niN as the up and down 
spin electron numbers respectively. That is, we consider a 
quantum state with conserved densities of up and down- 
spin electrons. Because of the conserved electron densi- 
ties, the single-site RDM of site i has a diagonal struc- 
ture, using the site basis states as unoccupied, doubly- 
occupied, up-spin-only occupied, down-spin-only occu- 
pied states, it is given as 



(3) 



In the above d is the probability of double occupancy, 
which related to the total number of doubly-occupied 
sites in the state as Nd = D =< ^ n.iin.ii). 

From the structure of the single-site RDM, we can infer 
a hierarchy of entanglement in many-electron states. The 
spin-only states with evert site being occupied by 

an electron with either up or down spin, correspond to 
the half-filled case of strongly-correlated electron states. 
Here, we have d = 0, n = 1, i.e. no double occupancy 
and no holes, effectively reducing the local Hilbert space 
to two, and the single-site RDM has only two eigenval- 
ues. The maximum entanglement occurs when both these 
eigenvalues are same, i.e. = = 1/2. Thus we have 
from Eq.3, 
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It should be pointed out that there will be many spin 
states with this entanglement e — 2/3. In fact, all many- 
electron states with the total spin zero, no double occu- 
pancy, and no holes will satisfy this. 

The second category comprises of states {|'02)} with 
no doubly-occupied sites but a nonzero fraction of holes, 
d — 0,n < 1. Since the site Hilbert space now is three 
dimensional, we have three nonzero eigenvalues of the 
single-site RDM. These states correspond to low-lying 
states of a large-U Hubbard model, or very strongly- 
correlated electron states with less than half filling elec- 
tron density. Here, the maximum entanglement occurs 
for 1 — n = ri| = ri| = 1/3, with the eigenvalues of RDM 
being equal to 1/3, we have 
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Thus existence of holes increases the entanglement. 

Finally, let us consider states with doubly-occupied 
sites as well. Maximum entanaglement occurs if all the 
four eigenvalues of the RDM shown above are all equal, 
i.e. d = 1/4, n-f = ni = 1/2. This corresponds to the 



half-filling case, except that doubly-occupied sites are al- 
lowed. Since the local Hilbert space dimension has been 
maximized to four, this gives the maximum entangle- 
ment, a global maximum. 



c(V') = £{d 



1 ^ 



1. 



(6) 



4' ' '2' 
In general there will be many such states with e = 1, 
since the only criteria we used are the conserved electron 
densities and the spatial uniformity. However, it will be 
useful to analyze various physical states such as metallic 
states, strongly-correlated states, superconducting states 
etc. Below we will examined these qualitatively different 
states. 

Let us consider the Nagaoka state[l3|, with d = 
0,Ne — N — 1, which is the ground state of C/ 00 
Hubbard model. This corresponds to a multiplct with 
the maximal total spin S = Sm0.x = N — 1/2. However, 
all the states in the ground state multiplet need not have 
the same entanglement. Consider the ground state with 
5^ = iV-l-Z, i.e. = N-l-l,Ni = l. The eigenval- 
ues of the single-site RDM are 1/iV, 0,1- (1 + 1)/N,l/N, 
giving an entanglement £(5*^ = N—l — l) — ^ ^^^^ 
It is easy to see that the maximum entanglement occurs 
for the state with I = N — 1/2, we have 
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This illustrates that the presence of a hole increases the 
entanglement from 2/3, which is the maximum value for 
the case of half filling with infinite Hubbard interaction. 

Let us now examine the behavior of the entanglement 
as a function of d, the average density of the doubly- 
occupied sites. A state with an optimal d is the uncorre- 
lated metallic state or Fermi ground state, which is con- 
structed as a direct product of single-particle momentum 

states, |F) = O '^It'^fcil^^' ^^'^^ single-particle k < kp 
state is occupied by two electrons, and the electron den- 
sity determines the Fermi wave vector kp- In the site 
basis, the state does not have a direct product structure, 
thus exhibits entanglement. It is easy to see that in the 
Fermi ground state, the double occupancy is given by 
d = n^/4, and e = 1 for the case of half filling n — 1. 
A popular state that incorporates a double occupancy as 
the on-site correlation strength is varied is the Gutzwiller 
statefll], given as 
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where A/" is a normalization factor, g determines the am- 
plitude for a state with double site occupancy. It is clear 
that g — 1 corresponds to the uncorrelated metallic state, 
where as the other extreme g = corresponds to no dou- 
ble occupancy or strongest spin correlations, or to J7 = cx3 
Hubbard model ground states. The double occupancy in 
the above state is related to the normalization factor, 

1 91og7V 



d{g) 



2N a log 5 



(9) 



3 




0.95 - 
0.9 - 

0.85 
0.8 

0.75 







n — 1, Hujd/Ef = 0.1 

n = l,hujD/EF = 0.2 

n = 0.5, Hujd/Ef = 0.1 

n = 0.5,hujD/EF = 0.2 

I I I I 

0.2 0.4 0.6 0.8 

Aq 



FIG. 1: The entanglement in the Gutzwiller state as a func- 
tion of g the projection amplitude for different electron den- 
sities. The metallic state is at g = 1 and the strongly- 
correlated state ai g — 0. For small filling fraction the 
strongly-correlated state has a slightly higher entanglement 
than the uncorrelated state, which itself is quite small com- 
pared to the global maximum of the matallic state at n = 1. 



FIG. 2: The entanglement in the superconducting state as a 
function of the energy gap Aq, for a few values of the electron 
density and the Debye frequency. The global entanglement 
measure decreases for a nonzero gap, indicating that the nor- 
mal state at Aq = with an optimal double occupancy max- 
imizes the entanglement. 



which is not easy to calculate in general. In one dimen- 
sion, it has been calculated by Metzner and Vollhardt [isj . 
we have 
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Using this, the entanglement is easily calculated, and is 
shown in Fig.l as a function of the projection parameter 
g. It is seen that optimal entanglement occurs for the 
uncorrelated metallic state at 5 = 1 at half filling, and for 
small electron densities the entanglement at near g ^ 
can be greater than the corresponding metallic case near 
g ^ 1, though, however, the actual amount is quite small 
compared to the global maximum at n = !,.'?= 1. Thus, 
the tendency of inhibiting double occupancy, from the 
optimal metallic value of d = n'^/i , decreases in general 
the site mixedness or the global entanglement. 

We now examine situations where it is energetically fa- 
vorable to promote double occupancy, over and above the 
optimal metallic value. Superconducting states fit into 
this category, as the underlying attractive interaction 
mechanism tends to promote superconducting Cooper 
pairs of opposite spins and in a zero momentum state, 
which is equivalent to a promoting tendency for site dou- 
ble occupancy. Starting with the vacuum state |0), the 
superconducting state can be constructed as 



Here, c[.| creates an up-spin electron in a single-particle 
plane-wave state with energy Ei^ — ti^k'^/2m. is the 
amplitude for a Cooper pair to form, Uk is fixed through 



the normalization lufel -|- \vk\ 



1 for every k. The 



Cooper pair probability is chosen as[l| 
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An attractive interaction in a narrow energy range near 
the Fermi energy Ep, would imply that the gap is nonzero 
only in that energy range, « Aq, |-Efc — -E_f| ~ fv-oo, 
where, uod '^s the Debye frequency or a suitable energy 
scale in general superconductors. 

The signle-site RDM is different from the displayed one 
in Eq.4, since the number of particles is not a conserved 
in this state, and is given by 
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In the above the off-diagonal matrix element is related 
to the probability of pairing on a given site, and can be 
written in terms of the Cooper pair amplitudes as, C, = 
(ci|Ci|) = ^UfcW-fc. It is straightforward to calculate, 
we have 
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The double occupancy in this state is given by d = 



(11) \ , which is larger than the optimal metallic value 
due to the pairing tendency as we argued earlier. Now, it 
is easy to calculate he global entanglement measure in the 
state, which is shown in Fig. 2 as a function of the super- 
conducting order parameter for a few values of Hlud/Ep 
and the electron density n. It can be seen from the plot, 
that for nonzero (which means the double occupancy 
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FIG. 3: The on-site concurrence between up and down spins 
in a superconductor for different electron densities and Debye 
frequencies. 

is more than the optimal metaUic value) the entangle- 
ment decreases, though only slightly for the half-filled 
case, about 15 per cent for smaller electron filling. 

The promotion of double occupancy by explicit intro- 
duction of off-diagonal correlations can introduce an en- 
tanglement structure of different variety, viz. a single-site 
concurrence, similar to the two-site concurrence in spin- 
only states Two-site concurrence in spin systems 
is a measure of the entanglement between two sites in 
a many-spin state, viz. the bipartite entanglement be- 
tween two spatial parts labeled by the two sites. Here, 
in the superconducting state, the one-site RDM can be 
viewed as a bipartite density matrix of two distinct spin 
parts. Let us rewrite the basis states of the RDM as 



|0 >= |0 >T |0 >x,| Ti>= |1 >T |i >i>l T>= |i >T 
|0 >i7 1 i>= |0 >| |1 >x- Now the basis states have the 
bipartite direct-product structure of up-spin states and 
down-spin states. Naturally, one can ask how much en- 
tanglement is there between the up spin and the down 
spin (which are used to label the two distinct parts of the 
bipartite Hilbert space). Using the Wootters' formula for 
he concurrence of the four-dimensional RDM, we have the 
single-site concurrence in the superconducting state as 

CTa=2Max(C-|--d|,0). (15) 

Thus, there is possibility of a nonzero on-site concur- 
rence, depending on the pair order parameter, or the en- 
ergy gap of the superconducting state. For example, we 
have C|,x ^ 0, if C(l + C) > f (1 - f)- We have plotted 
the single-site concurrence calculated above in Fig. 3 as 
a function of Aq /hujjj, for a few representative values of 
the fermi energy and the electron filling. 

In conclusion, we have examined the average impurity 
of the single-site reduced density matrix as a global en- 
tanglement measure in many-electron systems. The max- 
imum global entanglement occurs for the half filled uncor- 
related metallic state, which optimizes densities of up and 
down electrons, the holes and the doubly-occupied sites, 
"-T ~ "^i ~ i^hoie = d — 1/4. On-site density-density 
repulsive interactions, or strong spin correlations, tend 
to inhibit double occupancy, thus decrease the entangle- 
ment, as illustrated by Gutzwiller projection states. Su- 
perconductor states promote double occupancy, over and 
above the optimal value of the metallic states, but still 
the entanglement decreases. We have shown that there 
is possibility of a nonzero on-site concurrence in the su- 
perconducting states depending on the energy gap. 
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